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Abstract. We present 3D, gasdynamic simulations of jet/cloud collisions, with the purpose of modelling the 
HH 270/110 system. From the models, we obtain predictions of Ha and H2 1-0 s(l) emission line maps, which 
qualitatively reproduce some of the main features of the corresponding observations of HH 110. We find that the 
model that better reproduces the observed structures corresponds to a jet that was deflected at the surface of 
the cloud ~ 1000 yr ago, but is now boring a tunnel directly into the cloud. This model removes the apparent 
contradiction between the jet /cloud collision model and the lack of detection of molecular emission in the crossing 
region of the HH 270 and HH 110 axes. 
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1. Introduction 

The Herbig-Haro (HH) jet HH 110 (discovered by 
Reipurth & Bally 1986) is the best observed example of 
a possible HH jet/dense cloud collision. Reipurth et al. 
(1996) have interpreted the rather unique, collimated but 
quite chaotic structure of HH 110 as the result of a de- 
flection of the faint HH 270 jet through a collision with a 
dense cloud. 

The evidence presented by Reipurth et al. (1996) for 
this interpretation can be summarized as follows: 

— no stellar source has been detected aligned with the 
HH 110 jet, 

— the HH 270 jet (ejected from a detected IR and radio 
source, see Rodriguez et al. 1998) points towards the 
beginning of the HH 110 jet, 

— the proper motions of HH 270 and HH 110 have an 
approximately 2 to 1 ratio, which is completely con- 
sistent with the « 60° deflection angle defined by the 
locci of the two jets (the flow approximately lying on 
the plane of the sky) . 
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This last statement can be understood as follows. When 
a radiative jet hits the surface of a dense cloud at an in- 
cidence angle <fi (between the incident jet axis and the 
cloud surface), the normal component of the jet velocity 
is stopped in a radiative shock, and the jet continues to 
flow parallel to the surface (conserving the component of 
the incident jet velocity parallel to the surface). Therefore, 
the velocity Vdef of the deflected jet is approximately equal 
to the projection of the incident jet velocity Vi nc along the 
surface of the dense cloud. Therefore, v^ef ~ Vi nc cos <f>. 
One can straightforwardly see that the proper motions 
and deflection angle defined by HH 270 and HH 110 (see 
above) do satisfy this condition. 

This result led Raga & Canto (1995) to study the dy- 
namics of the collision of a radiative, HH jet with the sur- 
face of a dense cloud. These authors presented an analytic 
model and plane, 2D simulations of the early stages of such 
an interaction, and found that the general characteristics 
of the HH 270/110 system could be explained in terms of 
such a model. The main problem found with the models is 
that in a rather short timescale, the incident jet starts to 
perforate the obstacle, and the deflected jet beam is then 
pinched off. In order to obtain a long enough timescale for 
the production of the deflected jet, it is necessary to have a 
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Fig. 1. Time sequence of the density stratifications obtained from model A. The density stratifications on the y = 
plane (which includes the outflow axis and the centre of the spherical cloud) are shown for different integration times 
(as indicated at the bottom left of each plot). The densities are depicted with a logarithmic greyscale, with the values 
given (in g cm -3 ) by the bar on the top left of the figure. The x (horizontal) and z (vertical) axes are labeled in cm. 



very high cloud-to-jet density ratio. Raga & Canto (f 995) 
suggested that this problem might be overcome if the in- 
cident jet did not have a completely straight jet beam, so 
that the impact point would roam over the surface of the 
dense cloud. 

The regime in which the jet has punched a hole 
through a cloud was described by Canto & Raga (1996) 
and Raga & Canto (1996). If the cloud is stratified, the 
path of the jet through the cloud is curved, though the 



curvature is important only if the radius of the cloud is 
comparable to the jet radius. 3D gasdynamic simulations 
of the penetration of a jet into and through a dense, strati- 
fied cloud were carried out by de Gouveia Dal Pino (1999). 

Finally, Hurka et al. (1999) have studied the bending 
of the beam of a 3D MHD, non-radiative jet by a mag- 
netic field with a strong gradient (as would be found at 
the surface of a dense cloud). These authors show that 
this effect would help to increase the timescale over which 
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Fig. 2. Temperature stratification (top) and adaptive grid 
structure (bottom) on the y — plane obtained from 
model A for at = 2600 yr integration time. The tempera- 
ture stratification is depicted with a logarithmic greyscale 
with the values given (in K) by the bar on the top of the 
figure. In the bottom plot, two thick lines separating the 
jet, cloud and environmental material are shown (see the 
text). These lines show two values of the passive scalar : 
ip = (outer contour) and ip = 1.5 (inner contour). 



Fig. 3. Constant density surface (corresponding to a n = 
20 cm -3 number density of atomic nuclei) from model A 
for a, t — 2600 yr integration time. The two graphs show 
the surface as seen from two different directions. 



— that the jet is ejected with a constant direction and 
velocity 

— that it is produced with a precessing outflow direction 
and a sinusoidally varying velocity. 



the jet/cloud surface interaction takes place, before the 
deflected jet is pinched off. 

In the present paper, we discuss 3D gasdynamic simu- 
lations of the interaction of a radiative jet with the surface 
of a dense cloud. We show the results from two simulations 
with different assumptions for the incident jet : 



Through a comparison of these two simulations, we can 
evaluate the effect of a "roving" impact point on the pro- 
duction of the deflected jet. 

Our simulations include a treatment of the dissociation 
and ionization of the gas. Therefore, we can use the results 
to obtain predictions of atomic and molecular lines, which 
we directly compare with previously published images of 
HH 110. 



4 



Raga et al.: Jet/cloud collision simulations 



In particular, we compute predicted maps in the H2 
1-0 s(l) line. This is of interest because the morphology of 
HH 110 in this IR line is quite strikingly different from its 
morphology in atomic/ionic lines. Davis et al. (1994) and 
Noriega-Crespo et al. (1996) found that the H 2 emission 
is much better collimated, and lies along one of the edges 
of the HH 110 jet beam. This led Noriega-Crespo et al. 
(1996) to present a simple model of the molecular emission 
as coming from material from the dense cloud which is 
entrained by the jet as it brushes past the cloud surface. 
Our present simulations allow us to make a more definite 
assesment of whether or not such a mechanism actually 
succeeds in explaining the molecular emission of HH 110. 

We should point out that Choi (2001) present HCO+ 
emission maps, in which they detect emission in the 
HH 270/110 region, but not around the "point of impact" 
in which the "incident" HH 270 jet is presumably redi- 
rected into the "deflected" HH 110 jet. This result leads 
them to suggest that HH 110 might actually not be the 
result of a jet/cloud collision, but that it could instead be 
a "straight" jet ejected from a low luminosity, undetected 
stellar source which is presumably more or less aligned 
with the direction of the HH 110 flow. In the conlcusions, 
we discuss the possible ways of reconciling the jet/cloud 
interaction model with the observations of Choi (2001) 
which are suggested by our 3D gasdynamic simulations. 

2. The numerical simulations 

2.1. General features 

We have carried out 3D gasdynamic simulations of 
jet/cloud interactions with the yguazu-a adaptive grid 
code. This code integrates the 3D (or 2D) gasdynamic 
equations together with a set of continuity equations for 
atomic/ionic or chemical species. The details of the gas- 
dynamic and the adaptive grid algorithms have been pre- 
sented by Raga et al. (2000), and tests of the code are 
given by Sobral et al. (2000) and Raga et al. (2001). 

For the present simulations, we have used the follow- 
ing set of species : H2, HI and II, C II, III and IV, and 
O I, II and III (with abundances by number relative to hy- 
drogen : y c — 6.6 x 10~ 4 and y = 3.3 x 10~ 4 ). For the 
atomic/ionic reactions, we have included the collisional 
ionization (from Cox 1970), radiative+dielectronic recom- 
bination (from Aldrovandi & Pcquignot 1973, 1976) and 
O/H charge exchange processes. The H 2 dissociation and 
cooling has been included in the same way as in Raga et 
al. (1995). The cooling associated with the atomic/ionic 
species has been included as described in appendix A. 

We have computed two jet/cloud interaction models, 
which share the following characteristics. In both mod- 
els, an initially atomic jet (except for C, which is singly 
ionized) of number density rij = 50 cm -3 and temper- 
ature Tj = 1000 K interacts with a spherical, homoge- 
neous molecular cloud (with all H in the form of H 2 ) 
of number density n c = 5000 cm~ 3 and T c = 1 K. 
The cloud is surrounded by a homogeneous, neutral en- 



vironment of density d env = 10 cm -3 and temperature 
Tenv = 1000 K. Therefore, the jet to cloud (mass) density 
ratio is Pj/p c = 1/100. 

In both simulations, the cloud has a r c = 4 x 10 17 cm 
radius, and the jet has an initial, top-hat cross section 
of radius rj = 1.5 x 10 16 cm. The jet is injected at 
(x, y, z) = (0, 0, 0) with the outflow axis in the z-direction. 
Free outflow conditions are applied on all of the outer 
boundaries of the computational domain, except for the 
z = plane, on which a reflection condition is imposed 
outside of the injection jet cross section. 

The computations are carried out on a 5-level, bi- 
nary adaptive grid (the two coarsest levels being defined 
over the full computational domain, see Raga et al. 2000) 
with a maximum resolution (along the three axes) of 
3.91 x 10 15 cm. The highest resolution grid is only allowed 
in the regions occupied either by jet or by cloud material 
(which are traced by advecting a passive scalar), so that 
the maximum resolution allowed in the environmental ma- 
terial is only of 7.81 x 10 15 cm. 

We have then computed two models, one with a jet 
with time-independent injection conditions (model A), 
and one with a precessing, variable ejection velocity jet 
(model B). These two models are discussed in the follow- 
ing two subsections. 

2.2. Jet with time-independent injection (model A) 

In this simulation, the jet is injected parallel to the z- 
axis, with a constant Vj—300 km s _1 injection velocity. 
The computational domain extends from —4.5 x 10 17 < 
x < 0.5 x 10 17 cm, -2.5 x 10 17 < y < 2.5 x 10 17 cm and 
< z < 5 x 10 17 cm. The centre of the spherical cloud 
is placed at (x c , y c , z c ) — (3.5, 0, 3) x 10 17 cm, so that the 
jet has a glancing collision with the surface of the cloud. 

Figure 1 shows a time series (spanning an integration 
time of t — 2400 yr) of the density stratifications on the 
y = plane (this plane includes the outflow axis and the 
centre of the dense cloud). In this time series one sees the 
incident jet beam (injected at the origin, and travelling 
along the z-axis) impinging on the surface of the dense 
cloud, and being deflected onto a direction towards the top 
left of the xz-cuts. At t = 1200 yr, the jet has already dug 
a hole into the cloud (this hole becoming deeper at later 
integration times), and the deflected jet beam basically 
becomes cut off at its base. At this time, the jet/cloud 
impact point lies within, rather than at the surface of the 
cloud. However, the material deflected by the cloud surface 
at earlier times continues to travel away from the cloud, 
leaving a complex "wake" , joining it to the point at which 
the incident jet impacted the cloud surface. 

In order to illustrate the configurations adopted by the 
adaptive grid, Figure 2 shows the temperature stratifica- 
tion and the adaptive grid structure on the y = plane 
obtained for t = 2400 yr. It is clear from this figure that 
the higher resolution is not allowed on the regions occu- 
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Fig. 4. Time sequence of the density stratifications obtained from model B. The density stratifications on the y = 
plane (which includes the outflow axis and the centre of the spherical cloud) are shown for different integration times 
(as indicated at the bottom left of each plot). The densities are depicted with a logarithmic greyscale, with the values 
given (in g cm -3 ) by the bar on the top left of the figure. The x (horizontal) and z (vertical) axes are labeled in cm. 



pied by environmental gas, so that the leading bow shock 
is only resolved at the second highest resolution level. 

Figure 2 also shows the following. We have integrated 
an advection equation for a passive scalar ip. This scalar 
has been given a value of ip = 1 for the jet, 2 for the 
dense cloud, and —1 for the surrounding environment. In 
the plot showing the adaptive grid, we have also drawn 
two contours on the stratification of the passive scalar 
corresponding to values ip — (outer contour) and ip = 1.5 
(inner contour). The region in between the two contours 



corresponds to the jet material (which has ip = 1). From 
Figure 2 it is then clear that the jet material occupies 
the injection region and the hole in the cloud, as well as a 
"plug" of material (at (x,y) « (-2.5,4.5) x 10 17 cm) with 
wings which extend towards the jet/cloud impact point. 
The region in between the wings is filled in by a tongue 
of cloud material which has been swept into the deflected 
jet flow. We find that this entrained cloud material has 
interesting observational properties, which are described 
in section 3. 



6 



Raga et al.: Jet/cloud collision simulations 



The structure of this dense tongue is more clearly 
seen in Figure 3, which shows a constant density 3D sur- 
face (corresponding to a n — 20 cm~ 3 number density 
of atomic nuclei) obtained for a t — 2600 yr integration 
time. This hgure shows the jet penetrating into the cloud, 
part of the bow shock, the denser region of the deflected 
jet beam and the entrained molecular gas material (which 
forms a structure which surrounds the incident jet, and 
has an elongation in the direction of the deflected jet 
beam). 

2.3. P 'recessing, variable ejection velocity jet (model B) 

We have computed a second simulation, in which the ejec- 
tion direction precessses around the jet axis. The pre- 
cession cone has an a = 5° half-opening angle, and 
a t p — 400 yr precession period. Also, the jet is in- 
jected with a constant density (nj = 50 cm~ 3 , see §2.1), 
but with a sinusoidally varying ejection velocity Vj{t) = 
(300 + 80sin27rf/T„) km s _1 , with a t v = 200 yr period. 

For this simulation, we choose a computational domain 
with the same extent as the one of model A along the x 
and y axis (but with a diferent centering : —3.5 x 10 17 < 
x <1.5x 10 17 cm and -2.5 x 10 17 < y < 2.5 x 10 17 cm) 
but with a larger, < z < 10 18 cm extent along the axis. 
The center of the spherical cloud is placed at (x c , y c ,z c ) = 
(3.5, 0, 6) x 10 17 cm. We have chosen to have a larger dis- 
tance from the point of injection to the jet/cloud colli- 
sion region in order to allow the internal working surfaces 
of the jet (which result from the ejection velocity time- 
variability) to form before colliding with the dense cloud. 

Figure 4 shows a time sequence of the density strati- 
fications obtained on the y = plane (which includes the 
precession axis and the centre of the cloud). In this figure, 
one sees the internal working surfaces which first form, and 
then impact the surface of the dense cloud. Because of the 
precession in the ejection direction, the successive working 
surfaces impact the cloud at different points. Through a 
comparison with Figure 1, it is clear that this effect in- 
creases the time that the jet takes to dig a hole into the 
cloud (and therefore pinching off the deflected jet beam). 
In fact, for the t — 3000 yr time-integration shown in 
Figure 4, the depth of the hole is still smaller than the 
diameter of the impinging jet (for a similar time frame 
obtained from model A, the depth of the hole is of ap- 
proximately two jet diameters, see Figure 1). Similar re- 
sults are deduced by analyzing different y = const, cuts 
through the 3D density stratification. 

Figure 5 shows the temperature stratification, and the 
grid structure on the y = plane obtained for a t = 
2500 yr time-integration (equivalent results for model A 
are shown in Figure 2). On the graph with the adaptive 
grid, we again show the contours that separate the jet, 
cloud and environmental material (for model A, see Figure 
2 and the discussion at the end of §2.2). It is clear that 
model B has a more complex structure than model A, 
showing a number of condensations of jet material in the 
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Fig. 5. Temperature stratification (left) and adaptive grid 
structure (right) on the y = plane obtained from model 
A for at — 2500 yr integration time. The temperature 
stratification is depicted with a logarithmic greyscale with 
the values given (in K) by the bar on the top of the figure. 
In the right hand side plot, two thick lines separating the 
jet, cloud and environmental material are shown (see the 
text). 

"deflected flow" region, which correspond to the different 
working surfaces that have been deflected on collision with 
the cloud surface. Interestingly, the center of the deflected 
flow region is filled with material from the molecular cloud, 
as is also the case for model A (see Figure 2). 

In the following section, we present predictions of emis- 
sion line intensity maps carried out from the results of 
models A and B. These predictions can then be compared 
directly with the available observations of the HH 270/110 
system. 

3. Predicted intensity maps 

3.1. General considerations 

In order to compare the jet/cloud interaction models de- 
scribed in §2 with the published images of the HH 270/110 
system, we have obtained predicted emission line maps 
from the models. From the computed temperature, den- 
sity, electron density and H ionization and molecular frac- 
tions, we have computed the Ha and H 2 TO s(l) (2.12^m) 
emission coefficients. 

For the Ha emission line coefficient, we have included 
the recombination cascade and the n = 1 — > 3 collisional 
excitation (using the corresponding excitation coefficient 
of Giovanardi & Palla 1989). The H 2 1-0 s(l) coefficient 
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Fig. 6. Ha maps predicted from model A, corresponding 
to the integration times given on the bottom left of each 
frame. The maps are depicted with the greyscale given 
by the bar at the top of the graph (which gives the in- 
tensity values in erg cm~ 2 s _1 sterad -1 ). The two top 
frames have been computed including the dust extinction 
of the dense cloud (see §3.2), and the bottom frame shows 
a map computed without considering this extinction. The 
proper motions computed from the positions of three in- 
tensity maxima (as measured in the t = 2000 and 2600 yr 
maps) are shown in the top frame. The x (horizontal) and 
z (vertical) axes arc labeled in cm. 




Fig. 7. Ha and H2 1-0 s(l) intensity maps computed from 
model A for a t = 2600 yr integration time (see §3.2). The 
left plot shows the Ha map in factor of 2 contours and the 
H2 map in the greyscale described by the bar at the top 
of the plot. The right plot shows the H 2 map in factor of 2 
contours and the Ha map in gresycale (corresponding to 
the bar at the top of the plot). The greyscales (given in 
erg cm~ 2 s _1 stcrad -1 by the corresponding bars) and the 
contours of the intensity maps of a given line correspond 
to the same range of intensities. The x (horizontal) and z 
(vertical) axes are labeled in cm. 



was computed by solving the level population equations 
according to the prescription of Draine et al. (1983, cor- 
rected according to Flower et al. 1986). This is by no 
means the more accurate calculation of H2 level popu- 
lations available (see, e. g., Flower & Pineau des Forets 
1999), but it is appropriate given the limited accuracy of 
our rather low resolution numerical simulations. 

We have then computed intensity maps by integrating 
the emission coefficients along lines of sight. We have as- 
sumed that the y = plane (which includes the axis of the 
incident flow and the centre of the cloud) coincides with 
the plane of the sky. This is probably a reasonable approx- 
imation to the orientation of the HH 270/110 flow, since 
it is known that both HH 270 and HH 110 approximately 
lie on the plane of the sky (Reipurth et al. 1986). 

3.2. Intensity maps predicted from model A 

In Figure 6, we show the Ha emission line maps predicted 
from model A (see §2.2). From the bottom plot, it is com- 
pletely clear that the emission is dominated by the region 
in which the incident jet beam is digging into the cloud. 
As the optical emission maps of HH 110 do not show this 
emitting region, we have to conclude that it has to be 
absorbed by the dust in the dense cloud (at least, if we 
believe in the jet/cloud interaction model for this object). 

In order to illustrate the effect on the Ha maps of such 
an extinction, we have computed emission maps including 
this effect (top two maps of figure 6). We have computed 
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the maps assuming that the dust extinction coefficient is 
Kd = 10 _20 (n^/cm _3 ) cm -1 , giving a Td — 40 optical 
depth through the radius of the cloud. 

We should note that if we use the standard, n d = 
10~ 21 (n ff /cm~ 3 ) cm -1 dust absorption coefficient, our 
cloud would have a Td = 4 central optical depth. The 
region on the edge of the cloud into which the jet is pen- 
etrating would then have a low optical depth, producing 
little extinction of the emission from the impact region. 
This, however, is not a major problem given the fact that 
the dense cloud present in the HH 270/110 region is by no 
means either spherical or homogeneous, and could easily 
produce a large extinction towards the current jet/cloud 
impact region (which would lie within the cloud). 

The Ha maps obtained from model A (for integra- 
tion times t = 2000 and 2600 yr) computed with the dust 
extinction as described above are shown in the two top 
frames of Figure 6. It is clear that these maps do present 
a qualitative similarity to HH 110 (see, e. g., Rcipurth et 
al. 1996). 

In agreement with the observations, the emitting re- 
gion starts with a bright rim (in contact with the sur- 
face of the cloud, see figure 6) which points to a broader 
structure (with a complex structure of curved ridges) at 
larger distances from the impact region. Also, there is a 
faint emission halo extending parallel to the main emis- 
sion structure on the side directed away from the dense 
cloud. This is also in qualitative agreement with the Ha 
maps of HH 110 (Rcipurth et al. 1996). 

From the two time frames shown in Figure 6, we have 
computed proper motions for the three main intensity 
maxima seen in the maps. We have also computed proper 
motions for some of the local maxima in the region in 
which the emission is in contact with the cloud surface, 
but the resulting velocities lie between 3 and 8 km s _1 , 
and have not been plotted in Figure 6. 

The proper motions of the knots farther away from the 
cloud (shown in figure 6) have values of 15 to 45 km s _1 . 
These proper motions are substantially lower than the 
ones measured for the HH 110 knots, as Reipurth et al. 
(1996) have found values ranging from 35 to 150 km s" 1 . 

In Figure 7, we show a comparison between the Ha and 
the H 2 1-0 s(l) intensity maps obtained for t = 2600 yr. 
For computing the H2 map, we have considered an extinc- 
tion equal to 1/10 of the visual extinction. This results 
in only small optical depths towards the current jet/cloud 
impact region, so that the emission from this region is 
clearly visible. 

The H 2 emission in the base of the deflected jet re- 
gion is much more concentrated towards the surface of 
the cloud than the Ha emission. This result is in clear 
qualitative agreement with the morphology observed in 
HH 110 (Davis et al. 1994; Noriega-Crespo et al. 1996). 
At larger distances along the deflected jet flow, the H 2 
and Ha emission show spatially coincident condensations, 
again in agreement with the observations of HH 110. 



3.3. Intensity maps predicted from model B 

In Figure 8, we show the Ha intensity maps computed 
from model B for t = 2200 and 2400 yr integration times, 
and the H 2 1-0 s(l) map obtained for t = 2400 yr. Because 
of the fact that the jet/cloud impact region still lies on 
the surface of the cloud, the effect of the extinction due 
to the dust present in the cloud is not important, and 
has not been included. However, the extinction would be 
important for maps computed for different orientations of 
the jet/cloud structure with respect to the plane of the 
sky. 

As a result of the precession and ejection velocity time 
variability of the incident jet (see §2.3), the emission maps 
show more complex structures than the ones obtained 
from model A. In particular, one can clearly see the emis- 
sion from bow shocks around dense "bullets", which re- 
sult from the successive internal working surfaces present 
in the incident jet. Even though the morphology observed 
in the emission line maps of HH 110 is very complex (see, 
e. g., Reipurth et al. 1996), it does not appear to have such 
bow shock structures. Actually, the intensity maps pre- 
dicted from model A do resemble the structure of HH 110 
in a more convincing way. 

An interesting feature of model B is that the proper 
motions of the different intensity maxima (obtained by 
comparing the t = 2200 and 2400 yr Ha intensity maps, 
see figure 8) have values of close to 100 km s _1 . These 
velocities are in better agreement with the ones measured 
for HH 110 than the ones obtained from model A (see §3.2 
and Reipurth et al. 1996). 

4. Conclusions 

We have presented two jet/cloud collision 3D gasdy- 
namic simulations : one with an incident jet with time- 
independent injection conditions (model A), and a second 
one with a variable velocity, precessing incident jet (model 
B). A p c / Pj = 100 cloud to initial jet density ratio has 
been chosen for both models. 

Model A produces a deflection of the jet beam only 
for a ~ 500-1000 yr timescale, after which the incident 
jet starts digging a straight tunnel through the dense 
cloud. At later times, the deflected jet material contin- 
ues to travel away from the impact region, leaving behind 
a complex "wake" . 

Model B produces a broader jet/cloud impact region 
as a result of the jet precession. This effect results in a 
longer timescale for the duration of the jet deflection on 
the cloud surface (in fact, the jet deflection is still occuring 
at the end of our t = 3000 yr numerical simulation). 

Model B is more successful at reproducing the proper 
motions of HH 110, giving ~ 100 km s _1 velocities for 
the Ha intensity maxima along the deflected jet beam. 
Model A gives velocities of ~ 15-45 km s _1 , which are 
substantially lower than the proper motion velocities of 
HH 110 (see Reipurth et al. 1996). 
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Fig. 8. Intensity maps obtained from model B (see §3.3). The left and central frames show the Ha maps obtained 
for two integration times (given on the bottom left of each frame). The right hand side frame shows the H 2 TO s(l) 
map obtained for t = 2400 yr. The maps are depicted with the greyscales given (in erg cm" 2 s" 1 sterad" 1 ) by the 
bars on the top of each frame. The left plot shows the proper motions of several intensity maxima, computed from the 
t = 2200 and 2400 yr Ha intensity maps. The x (horizontal) and z (vertical) axes are labeled in cm. 



This difference between model A and model B is due 
to two effects. The first effect is that in model A, the jet 
is deflected only for a <~ 500-1000 yr timescale, and that 
this deflected jet material then slows down as it interacts 
with surrounding, environmental gas. In model B, succes- 
sive deflected "bullets" (i. e., internal working surfaces) 
travel into the low density region left behind by the pas- 
sage of the head of the deflected jet, and do not interact 
directly with the higher density environment. The second 
effect is that because of the precession of model B, some of 
the bullets have trajectories which are more tangential to 
the surface of the molecular cloud. These more tangential 
bullets are less deflected, and therefore preserve larger ve- 
locities than the ones that have a more normal incidence 
on the cloud surface (or than the deflected jet of model 
A). 

However, in most other counts, model A is more suc- 
cessful than model B at reproducing the observations of 
HH 110 : 

— the general qualitative appearance of the Ha maps is 
in better agreement with the HH 110 Ha images, 

— the features of the H 2 1-0 s(l) emission and their spa- 
tial relation to the Ha emission also are in good qual- 
itative agreement with HH 110, 

— the maps that better resemble HH 110 correspond to 
times (e. g., the t = 2400 yr frame of figures 1 and 6) 
in which the impact region is already immersed within 
the cloud. 



This last feature offers an interesting way of reconciling 
the jet/cloud collision model with the HCO + observations 
of Choi (2001). 

In these observations, no HCO + emission was detected 
in the region in which the axis of the "incident" HH 270 
jet crosses the axis of the "deflected" HH 110. Choi (2001) 
noted that this appeared to be in disagreement with a 
jet/cloud collision model for this system, as the cloud 
shock produced in the impact region should indeed pro- 
duce HCO + emission. 

The situation found in model A, however, could indeed 
be in agreement with the observations of Choi (2001). In 
this model, the impact region does not lie in the point 
in which the incident and deflected jet axes cross, but 
is instead located further along the axis of the incident 
jet. Interestingly, Choi (2001) does find substantial HCO + 
emission West of HH 110, approximately aligned with 
HH 270. 

As we have discussed in §3.2, the fact that the impact 
region is not observed optically can in principle be a result 
of the dust extinction in the dense cloud. Interestingly, 
some H 2 TO s(l) emission is apparently detected to the 
West of HH 110 (more or less aligned with HH 270, see 
Noriega-Crespo et al. 1996), which in principle might be 
associated with the impact region. 

To conclude, we note the two main features of our re- 
sults : 
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— we find that the jet/cloud interaction model does re- 
produce the Ha and H 2 1-0 s(l) emission observed in 
HH 110 in a qualitatively successful way, 

— our models show that the lack of a detected im- 
pact zone in the incident/deflected jet crossing regions 
(Choi 2001) is not a major problem for a jet/cloud 
collision model (as this region could presently be dis- 
placed further into the cloud). 

Clearly important questions remain about the more 
technical aspects of our simulations. In particular, the 
H2 emission from the deflected jet comes from molecular 
cloud material which has been entrained into the jet flow 
(see figures 2 and 5). Even though we find that our models 
produce H2 emission structures in agreement with the ob- 
servations of HH 110, the accuracy of our rather low reso- 
lution simulations in reproducing the cntrainment process 
(which gives rise to this emission) is somewhat question- 
able. Because of this, the present results have to be taken 
with some caution. 
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Appendix A: The atomic/ionic cooling rates 

For H I and II, C II and III and O I and II, we take the 
cooling rates (per atom and per electron) tabulated by 
Raga et al. (1997). In order to simplify the calculation 
of the cooling rate, we note that the cooling per atom 
or ion starts to deviate from the low density limit only 
for electron densities n e > 10 4 cm~ 3 . As such electron 
densities are generally not found in HH objects, one can 
safely consider only the low density limit. 

In this limit, the cooling per atom and per electron is 
only a function of the temperature, and one can then make 
simple parametric fits for the cooling due to each of the 
considered species. The following fits (giving the cooling 
rate in erg s" 1 cm 3 ) have errors smaller than 10% over 
the full 10 3 — > 10 6 K range (unless otherwise noted) : 

1. Collisional excitation of HI: 

login ( -^-) = -50+32.3(l-t)-1180(l-t°- 0001 ),(A.l) 
\n e n HI J 

with t = 1590K/T, 

2. Collisional ionisation of HI: 

LionMi = n e n HI q(T)xH , (A. 2) 

where q(T) is the collisional ionisation coefficient and xh 
the ionisation potential of H. 

3. Radiative recombination of HII: 

We have included the classical interpolation formula of 
Seaton (1959). 

4. Collisional excitation of O I (electrons): 

L = L\ + L 2 with 

login I Ll I = -23.95 + 1.23ii + 0.5ti 10 , (A.3) 
\n e noi J 

login f — ) = -21.05 + 1.2t 2 +1.2[max(t 2 ,0)] 5 ,(A.4) 

\n e n i J 

with h = l- 100 K/T and t 2 = 1 - 10 4 K/T. 

5. Collisional excitation of O I (neutrals) 

login! ] = 10.3t + i 8 -34.4, (A.5) 

\n HI n iJ 

with t = 10 K/T. This interpolation formula fits the com- 
puted cooling in the 10 — > 10 5 K temperature range. 
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6. Collisional excitation of O II 
L = L\ + L 2 with 



(—) 

\n e n n J 



log 10 — = 7.9h - 26.8 , (A.6) 



log 10 f— =1.9 20.5, (A.7) 

\n e nonJ N - 5 



with h = 1 - 2000 K/T and t 2 = 1 - 5 x 10 4 K/T. 



7. Collisional excitation of C II 

lo Sio ( Ll ) = 

- 23.65 + 1.2 [max(t - 2 ), 0] [1 - 5 -°- 25max( ^ 4 ' 0)1 , (A.8) 
with i = log 10 T [K]. 



8. Collisional excitation of C III 

log 10 f — — — ) = -20.8 + 3. 9(t-4)-1.37(t-4) 2 ,(A.9) 
\n e n C iii J 

with t = min(log 10 T [K], 5.4). 



9. Parametrized cooling for OIII 

We have replaced the real cooling due to collisional ex- 
citation of OIII with a function that resembles the coro- 
nal ionisation equilibrium for T > 10 5 K. In this way, the 
computed cooling has a transition to the coronal equilib- 
rium cooling when all O becomes OIII and all C becomes 
CIV (with which no cooling is computed). The adopted 
functional form is L = L\ + L 2 with 

log 10 ( ) = -21.4 - 3.5(t - 5.4) 2 , (A.10) 

\n e nyoinJ 

log 10 ( — ) = -21.7 - 0.7(t - 6.5) 2 , (A.11) 

\n e nyoinJ 

where t = log 10 T [K] and yoni is the OIII ionisation 
fraction. 



